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Мета роботи: придбання практичних навичок в чисельному 
інтегруванні звичайних диференційних рівнянь явними і 
неявними методами  Рунге-Кутта, дослідження впливу 
значення кроку обчислень на точність і збіжність рішення. 
Визначення можливості застосування засобів стандартних 
пакетів для отримання результатів. 


ПОРЯДОК ВИКОНАННЯ ЛАБОРАТОРНОЇ РОБОТИ 
1. Запрограмувати на мові пакету Маїпетайса рішення 
заданого диференційного рівняння явним методом Рунне - 
Кутта четвертого порядку (10.6) і виконати рішення при 
кількох значеннях кроку, поки рішення не почне розбігатися. 
Спробуйте з'ясувати, чи існує аналітичне рішення задачі. 
2. Порівняти отриманий максимально можливий крок Птах з 
значеннями, обчисленим за допомогою формули (10.17).. 
З. Запрограмувати на мові пакету Маїіһетаійса рішення 
заданого диференційного рівняння неявним методом Рунне - 
Кутта 4(5) і виконати рішення при максимальних значеннях 
кроку з пункту 1. 
4. Запрограмувати на мові пакету Маіһетайса рішення 
заданого диференційного рівняння вкладеним явним методом 
Рунне - Кутта четвертого порядку (10.8) і виконати рішення 
при максимально можливому кроці һтах ,знайденому в 
пункті 2. 
о. Користуючись стандартними операторами пакету 
Ма!етайса, знайти рішення заданого диференційного 
рівняння вкладеним явним методом Рунне - Кутта і порівняти 
покрокові похибки рішень, отриманих в пунктах 1, 315. 
6. Користуючись стандартними операторами пакету 
Ма!етайса, знайти рішення заданого диференційного 
рівняння неявним методом Рунне - Кутта Побудуйте 
відповідні графіки у часі. 
7. Скласти звіт з отриманих результатів і математичних 
формул використаних методів по кожному пункту завдання, 
давши оцінку порівняльної точності отриманих рішень 
різними методами. 


Варіант 24 


у'={(Ї, у) їо к Уо 


у21(2)/# 1 З 1 


Хід роботи 

1. Запрограмуємо метод Рунге-Кутта в Маіћетабса. 
Введемо межі відрізку інтегрування і, початкові умови уо, 
кількість кроків т, довжину кроку П та масив ( для 
зберігання результату розрахунків. 

не  хі:губьгодіві / Є: 

кі ]:= #([6, у]; 

когь різ #[+0.5+һ, у+0.5аһаКк1[һ]]; 

КЗ[һ | := [+ 0.5+һ, у+0.5+һғк2[Ь])];: 

ка[ћ | іс ЕЕ б, у+һ+Кк3[һ])];: 

0-1; 

са З; 

ОХ Е 

ш = 10; 

Һ = (їм - +0) / №; 

О = Аггау[о, (2, шо. 1), (0, 0)}]; 

аг0, 0] = +0; 

а(1, 0) = уо; 

е= 0: У= уб; 

ро {ду = пе (К1{һ) + 2+ К2(һ) + 2 + КЗ) + К4([Һ)) /6; 

= ++; у= у ду: о[0, 1) = 6; а(1, 1] = у, (1,1, ш)): 

Виведемо результати: 

Тар1еРогши [Тар1е[ (а10, 1], агї1, 1)), (1, 1, 10)), 

Тар1еНеадіпдз ~» (Айасошщабіс, ("С", "у"))) 


З 
6 
5 
5 
: 
5 
А 
5 

2 
1 


7 -. 1.62131 


8 з 1.83991 
9 = 2.12788 
10 | з 2.52184 


2. Обчислимо максимально можливий крок 


1(с ,у )-у"2аї1о6с(с)/6; 
1112 , у] «РІЄ(е, у), У] 
2 уІод!с) 


- 
С 


ЕіпаМахішом | (5112, у], 126 х«3), {{+, 1), (у, 1}}] 


132.4751, ‘с -2.71828, у з 44.1382)) 


Ћ = 2.7853 / 32.4751 


0.0857672 


Повторимо обчислення, але уже зһ = 0.085. 

ЕГЕ, У 1:- у"2аїо0саі БІ Е; 

кіІїь ]:= #[Е, у]; 

к2[һ | := Ғ[Е+0.5+һ, у+0.5+һаКк1[Ь]]; 

кЗ[А ] := #[Е+0.5+һ, у+0.5+һ+к2[һ]]; 

каїь різ Ғ[Е+Ь, уз пекзуіві): 

0-1; 

ба = 3; 

уФ-1; 

а = 23; 

ћ = 0.085; 

О = Аггау[0о, (2, ш. 1), (0, 0)]; 

аго, 0] = +0; 

0[1, 0] = уб; 

Е = 0; у = уб; 

рогду = Һ + (К1[Ь] +2+к2[] +2 + к3[Ь] + КА[Ь]) /6; 
Е= +0: узузФбу: піб, 11-55: сії, 11- т, Ті, 1; щі); 

Тарб1еРоги [Тар1е[ (аг0, 11, аг1, 11), (і, 1, 108), 
ТарієНеадіпдязв > (Апсошабіс, { "і", "у"))) 


беРогтє 


т 

1. 
1. 
1. 
1. 
1. 
1. 
22 
1. 
1. 
1. 


ы юу олњом№ ю 
ыыы ыыьыьыьыь ю ы 


Значення сильно відрізняються від попередніх. Можна 
побачити відмінність між отриманими з різним кроком 
розв'язками для однакових значень ї вже в третьому знаку 
після коми. 


З. Виконатиємо рішення при максимальних значеннях 
кроку з пункту 1. 


Ғ[Е ,а1с:-а?"Заїо0сібІ/Є; 
ћ = 0.085; 
ю = 10; 
О = Аггау а, т +1]; 
Т = Атгау[е, ваз 11: 
[0] =1; ұю] = 3; 
[01 = 21; 
ро [{+ [п] = [0] +һ+п; 
К1 := Ё[о[а+1], є[а+1]]; 
к2 ғ[ъ[а +1] -/2, а|по1)- (1/2) ++к1]; 
кЗ Е[+[2 +1] -Һ/2, а|по1) – (1/2) *+0к2]; 
КА := Е[[п0 +1], о[а +1] - бекз3і; 
а[0 +1] = о[п] +/6+ (к1+2+к2+2+к3 + к4)}, (п, 0, м}]; 
поь ГрезнЕгГен-е. в." Єр, п. 21-е. Ера. е геа. 7127, гат. фа, 0, аут: 


понт 


п=0 0[0]=1. 0[0]=1 
п=-1 0[1]=1.085 0[1]=1. 
п=2 0[2]=1.17 0[2]=1.00639 
п=3 02[3]=1.255 013]=1.01794 
п=4 0[4]=1.34 0[4]=1.03388 
п=5 0[5]=1.425 0[5]=1.05373 
паб 0[6]-1.51 0[6]=1.07719 
па? 0[7]=1.595 0[7]=1.1041 
п=8 0[8]=1.68 0[8]=1.13443 
п=9 —[9]=1.765 0[9]=1.16821 
п=10 С(101-1.85 ц(101-1.20555 
4. Розв'яжемо рівняння вкладеним явним методом Рунне - 
Кутта четвертого порядку і виконаємо рішення при 
максимально можливому кроці птах, знайденому в 
пункті 2. 
Для підвищення обчислювальної ефективності чисельного 
інтегрування диференціальних рівнянь і контролю покрокової 
похибки були розроблені методи, названі вкладеними 
методами Рунге-Кутта. Основна ідея полягає у виконанні 
кроку від точки Ё, до ї,,: двома методами порядку точності рі 
р + 1, для того щоб отримати оцінку похибки інтегрування 
методом порядку р. 
Найбільш розповсюдженим на сьогодні методом Фельберга є 
метод, який поєднує методи Рунге-Кутта четвертого і п'ятого 
порядків і має похибку апроксимації порядку т“ іт? відповідно. 


(6, У] :=у^20100[6) / Е: 
к2({( сац ) := 2[6 + 0.25 а гац, у + 0.25 а сац а к1([бац));: 
КЗ{ бац ) := 2(6+3/8а сали, у+3/32 а бапа к1 [ац] + 9/32 а сацак2(|сац)); 
К4 {баш ): 
2(6+ 12/13 а сац, 
у + Сац ж (1932 / 2197 ж К1({ бац] + 7200 / 2197 . К2 [ бац] + 7296 / 2197 а КЗ[ёаи])); 
К5{ бац ) := 
(Є + Саи, 
у + сап (439/216 а к1 [ац] - 8 а К2 [ бац] + 3680 / 518 + КЗ [бац] - 845 / 4104 + ка [ ёац])]; 
Ҝ6{ бац ) := 
2([# + 0.5 а ёли, 
У + 
сац ((-8/27) «Кі|сацп) - 2а к2({ сац) - 3544 / 2565 + КЗ{ сац) + 1859 / 4104 + к4 (сац) - 
211/40«Кк5|сацпі)): 
ш = 10: +0 = 1; у0 = 1; ба = 3; саа- 0.085; 
О = Атггау[а, {м +1, 2), 0): а(0, 0) = М(С0): 0(0, 1) = у0: 
Ер = Аггау [ер, {т +1, 2), 0]; 
ер(0, 0] = 10: ер[0, 1) = 0: 
3=0: = 0: У = уб: 
рої 
ау1 = 
Н ао (16 / 135 + к1 [бао] + 6656 / 12825 а КЗ [бап] + 28 561 / 56 430, к4 [ац] - 
9/50 + к5[+аџ] «2/55аКк6|саці)): 
ду = М[ ао + (25 / 216 а К1 (бао) + 1408 / 2565 + КЗ [бао] + 2197 / 4104 + К4[ ао) - 1/ 5аКк5(|саа))): 
ер(1, 1] = ду - 0у1; у = у + ау; о[і, 1) = у: С о М[Є о бао]; 
ер(і, 0) = 6: еёб[1, 0) = є: аїі, 0) = Є, (і, 1, в); 
М = Абз|М(Тарбіе(| (а(1, 0), а(1, 1), ер(1, 1)), (1, 0, 10)), 4)): 


ТарбіеРКогші|її, ТарієНеадіпдає > (Апасошабіс, { "і", "ша", "ерр"})}] 


с от ерр 

1 1.000 [2] 

2 1.085 1.00341 5.79164х 1078 
з 1-17 1.01319 0.0000549135 
4 11-295 1.02868 0.000129759 
5 1:34 1.0495 0.000219506 
6 1.425 1.07549 0.00032033 
7 1.51 1.10669 0.000431974 
8 1.595 1.1433 0.000556286 
9 1.68 1.18571 0.000696684 
10 1.765 1.23452 0.000858112 
11 1.85 1.29053 0.00104736 


5. Користуючись стандартними операторами пакету 
МаШетайса, знайдемо рішення заданого 
диференційного рівняння вкладеним явним методом 
Рунне - Кутта. 


зѕо1 = МО5оїуе| {у ' [1] == у[“] ^2+109[4] /ї, у[1] 22.1), У, (2,1, 3), 
Меёһоа з» "Ехр1ісі+ВопсдеКої+а"] 


у э ІпсегроїіатсіпдЕопсгіоп | Оотаіп: ((1., 3.) Пп 
Оџшіриё ѕсаіаг з 
УГЕ | с УГЄ) /. зо1 
ІпсегроїіасіпаКипссіоп || Оотаіп: (01., 3.) 1ге] 
Эиїриё ѕсаіаг 


Тао = 0.085; 
У = Тар1еЕогт [Тар1е[{1 + і + аъ, у[1 + і ғ саа) [[1]]}, (1, 1, 15)), 
Тар1іеНеааіпаз э (Апсошабіс, { "і", "у1"))) 


1 1.085 1.00334 
2 т 5 1.01248 
з 1.255 1.02648 
4 1.34 1.04474 
5 1.425 1.06691 
6 1.51 1.09279 
7 1.595 1.12232 
8 1.68 1.1555 
9 1.765 1.19246 
10 1.85 1.23339 
11 1.935 1.27855 
12 2.02 1.32833 
13 2.105 1.38312 
14 2.19 1.44352 
15 2.275 1.51017 


З отриманих результатів можна зробити висновки, що 


розв'язки, отримані явним методом 


Рунге-Кутта 


відрізняються від точного розв'язку на деяких ітераціях лише 


на 107. 


6. Користуючись стандартними операторами пакету 
Ма!етайса, знайти рішення заданого диференційного 
рівняння неявним методом Рунне - Кутта Побудуйте 


відповідні графіки у часі. 


801 = МО5о01уе|(у ' [+] - у[+] ^2 +109[+] / +, у[1] - 1), у, {+, 1, 3), 
Меёһћоа з "Ішр1ісі ВопсдеКоёёа"] 


[ Гу - ІпсегроїасіпаКкипссіоп 1 чинни Ма. 3) | } ) 
І: Оџшіриё зсаїаг 


УГЕ 1 с У[Є] /. ѕо1 


Оотаіп: {1., 3.} |]... 1 
чари Е 
Оџшіриё ѕсаіаг 


+аџ = 0.085; 
У = Тар1еЕогт [Тар1е[{1 + і + ап, у[1 + і ғ ао] [[1]]}, (1, 1, 15)}], 
Тар1іеНеааіпдѕ » (Апасошабіс, { "+", "у1"}}] 


Е у1 


Я 1.085 1.00334 
2 1-37 1.01248 
з 1.255 1.02648 
4 1.34 1.04474 
5 1.425 1.06692 
6 1.51 1.0928 
7 1.595 1.12232 
8 1.68 1.1555 
9 1.765 1.19246 
10 1.85 1.23339 
11 1.935 1.27856 
12 2.02 1.32832 
13 2.105 1.38313 
14 2.19 1.24352 
15 2.275 1.51017 


На деяких ітераціях останній розв'язок відрізняється від 
отриманих раніше приблизно на 107. 


Висновок 
Однокрокові методи Рунге-Кутта будують на основі 
скороченого ряду Тейлора. Кількість членів цього ряду р 


визначає порядок методу, а перший неврахований член ряду 
- локальну похибку обчислень О(ћ?*!) . Порядок методу р 
визначає залежність похибки від кроку обчислень. Чим вище 
порядок методу, тим точніше розв'язок, але й вище порядок 
обчислень. 

Недоліком явних методів Рунге-Кутта є те, що існує 
обмеження для максимального кроку обчислень, 
перевищення цього максимального кроку призводить до 
розбіжності розв'язку. Це робить явні методи непридатними 
для о розв'язання жорстких задач. В таких випадках 
використовуються неявні методи. 


